DESY 02-134 
hep-ph/0210161 
October 2002 



ISSN 0418-9833 



Heavy-quarkonium creation and annihilation with 

0{al In as) accuracy 

Bernd a. Kniehl, Alexander A. Penin* Matthias Steinhauser 

II. Institut fiir Theoretische Physik, Universitat Hamburg, 
Luruper Chaussee 149, 22761 Hamburg, Germany 

Vladimir A. Smirnov 

Nuclear Physics Institute, Moscow State University, 
119899 Moscow, Russia 



Abstract 

We calculate the 0{a^ Ina^) contributions to the heavy-quarkonium production 
and annihilation rates. Our result sheds new light on the structure of the high- 
order perturbative corrections and opens a new perspective for a high-precision 
theoretical analysis. We also determine the three-loop anomalous dimensions of the 
nonrelativistic vector and pseudoscalar currents. 
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The theoretical study of nonrelativistic heavy-quark-antiquark systems is among the 
earhest apphcations of perturbative quantum chromodynamics (QCD) [1] and has by 
now become a classical problem. Its applications to bottomonium [2] and top-antitop 
[3] physics entirely rely on the first principles of QCD. These systems allow for a model- 
independent perturbative treatment. Nonperturbative effects [4] are well under control for 
the top-antitop system and, at least within the sum-rule approach, also for bottomonium. 
This makes heavy-quark-antiquark systems an ideal laboratory to determine fundamental 
parameters of QCD, such as the strong-coupling constant and the heavy-quark masses 
ruq. The bottom-quark mass m(, is of particular interest, in view of current and future 
S-physics experiments. In the observables employed to extract the Cabibbo-Kobayashi- 
Maskawa matrix elements and to gain deeper insight in the nature of CP violation, m;, 
enters as a crucial input parameter [5]. Thus, precise knowledge of is essential for the 
interpretation of the experimental data. On the other hand, the top-quark mass rrit is one 
of the key parameters in the precision tests of the standard model of the electroweak inter- 
actions and in the search for new physics at a future e+e" linear collider. Furthermore, the 
study of tt threshold production should even allow us to probe Higgs-boson-induced effects 
[6]. Besides its phenomenological importance, the heavy- quarkonium system is also very 
interesting from the theoretical point of view because it possesses a highly sophisticated 
multiscale dynamics and its study demands the full power of the effective-field-theory 
approach. Equipped with reliable perturbative results and experimental data on heavy- 
quarkonium observables, one can test the effects and structure of the nonperturbative 
QCD vacuum. 

The binding energy of the bound state and the value of its wave function at the 
origin are among the characteristics of the heavy-quarkonium system that are of pri- 
mary phenomenological interest. The former determines the mass of the bound-state 
resonance, while the latter controls its production and annihilation rates. Recently, the 
heavy-quarkonium spectrum has been computed through 0{a^gmq) [7,8] including the 
third-order correction to the Coulomb approximation. On the other hand, as for the 
wave function at the origin, a complete result is so far only available through 0{a1) 
[9]. The 0(q;^) correction has turned out to be so sizeable that the feasibility of an 
accurate perturbative analysis was challenged [10], and it appears indispensable to gain 
full control over the next order. Only the double- logarithmic third-order correction, of 
0{al In^ as), is available so far [11]. In this Letter, we take the next step and calculate the 
single-logarithmic ^(a^lnas) correction. As a by-product of our analysis, we obtain the 
three-loop anomalous dimensions of the nonrelativistic vector and pseudoscalar currents, 
which constitute central ingredients for the renormalization-group improvement of the 
effective theory of nonrelativistic QCD (NRQCD) [12-14]. The main results are given by 
Eqs. (6), (7), and (10). As for the calculation, we follow the general approach of Ref. [7] 
(sec also Ref. [15]). It is based on the nonrelativistic effective-theory concept [16] in its 
potential-NRQCD (pNRQCD) incarnation [17] implemented with the threshold-expansion 
technique [18]. 

Let us focus on two examples of paramount phenomenological relevance: the Icptonic 
decays of the T(15') resonance and the threshold production of top quark-antiquark pairs 
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in e'^e~ annihilation. Both processes are essentially photon mediated and thus governed 
by the electromagnetic quark current = qji^q. Within the effective theory, has the 
following decomposition in terms of operators constructed from the nonrelativistic quark 
and antiquark two-component Pauli spinors ip and x [16]: 

j, = c,(/.)V'V,x + + ■ ■ ■ , (1) 

where /x is the renormalization scale, D are the space components of the gauge-covariant 
derivative involving the gluon fields, and the ellipsis stands for operators of higher mass 
dimension. The Wilson coefficients c^(/x) and dv{ii) may be evaluated as series in q;s(a*) 
and represent the contributions from the hard modes (where energy and three-momentum 
scale like ruq) that have been integrated out. They are computed in full QCD for on-shell 
on-thrcshold external (anti)quark fields and are logarithmic functions of ji/niq. Also 
integrating out the soft (energy and three-momentum scale like m^v, where v is the 
heavy-quark velocity) modes and the potential (energy scales like rngv"^, while three- 
momentum scales like rUgv) gluons yields the effective Hamiltonian of pNRQCD, which 
contains the potential (anti)quarks and the ultrasoft (energy and three-momentum scale 
like rUgV^) gluons as active particles. The dynamics of the nonrelativistic potential heavy- 
quark- antiquark pair in pNRQCD is governed by the corresponding effective Schrodinger 
equation and its multipole interactions with the ultrasoft gluons. The effective-theory 
expression for the partial decay width of T(IS') — > /+/^ reads [9] 



(J2 Q,2 

cKrrib) + ' cy{mb) {dy{mb) +3) + ... 



(2) 



with T^^ = 47rA^e<96a^|^f(0)|V(3m^) and pi = |V'i(0)|V|V'f (0)P, where N, = 3, is 
the fractional electric charge of quark q, a is Sommerfeld's fine-structure constant, ipi{x) 
is the ground-state wave function as computed in pNRQCD, and ipi{x) is the Coulomb 
solution, which incorporates the leading binding effects and about which the perturbative 
expansion of ipi{x) is constructed. For arbitrary principal quantum number n, we have 

^^(o)p = Cl^a^m^/ (87rn^), where Cf = {N^ - l)/(2A^c)- Here and in the following, 
Oisdi) is to be evaluated at the soft normalization scale fig = CFOis{lJ's)'m'q whenever its 
argument is omitted. Nonperturbative contributions to Eq. (2) are ignored. The leading 
one, due to the gluon condensate of the vacuum, may be found in Ref. [19]. It is quite 
sizeable and out of control for higher resonances. A reliable quantitative estimate of the 
nonperturbative contributions to Eq. (2) can only be obtained through lattice simulations. 
On the other hand, to keep the nonperturbative effects under control, one can employ 
nonrelativistic T sum rules [2] based on the global- duality concept. 

In the top-quark case, the nonperturbative effects arc negligible. However, the effect 
of the top-quark total decay width Tt has to be properly taken into account [3], as it 
is relatively large and smears out the Coulomb-like resonances below threshold. The 
NNLO-^ analysis of the cross section [10] shows that only the ground-state pole gives rise 

^In the effective-theory framework, one has two expansion parameters, Ug and v, and the corrections 
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to a prominent resonance. The value of the normahzed cross section R — a{e'^e~ — > 
ti)/(j{e~^e~ — > ii^jjr) at the resonance energy is dominated by the contribution from the 
would-be toponium ground-state, which is of the form 



Ri — Ri'^Pi 



(72^2 

{mt) + -^^v{mb) {dy{mb) + 3) + 



(3) 



with R\'~' = GnNcQjli/Ji {0)1"^ / (m^Ff). The contributions from the higher Coulomb-hke 
poles and the continuum are not included in Eq. (3), and we postpone the complete 
analysis to a future publication. It is understood that a appearing in F^^ and Ri'^ is to 
be evaluated at the mass scale of the respective resonance. 

Starting from c^ifJ^) is infrared (IR) divergent. This divergence arises in the 

process of scale separation and is canceled against the ultraviolet (UV) one of the effective- 
theory result for the wave function at the origin. In our approach, dimensional regular- 
ization with d — A — 2e space-time dimensions is used to handle the divergences, and the 
formal expressions derived from the Feynman rules of the effective theory are understood 
in the sense of the threshold expansion. This formulation of effective theory possesses two 
crucial virtues: the absence of additional regulator scales and the automatic matching of 
the contributions from different scales. For convenience, we subtract the IR and UV poles 
in e according to the modified minimal-subtraction (MS) prescription and set /i — niq, so 
that Cv{mq) is devoid of logarithms. The latter is known through 0{al) and reads [20] 



TT 



71 



+ 



897r^ 



^In2-^C(3) U7^C^ + 



151 
72" 
'23 797r2 



+ ..., 



36 



+ ^Mn2-ic(3))c|+(f -^)c^r^ 



(4) 



where as is renormalized in the MS scheme, Ca — Nc, Tp — 1/2, ni is the number of 
light-quark flavors, and C,{x) is Riemann's C, function with value ^(3) = 1.202057. . .. To 
the order considered, we have dy{mq) = 1. 
The corrections to 1^^1^(0)1^ read 



Pi 



+ 



as 
1 + — 

TT 



a.. 



4- — l/?o + ^ai 



TT 



cl 



are classified according to the total power of ag and v as leading order (LO), next-to-leading order (NLO), 
NNLO, N^LO, etc. 
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X TT^ HCpas) + f -— + 20C(3) + - j /?o' 



+ 4 



97r2 / 337^2 

+ — CaCf + 



137r^ 
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S{S + 1) C| 



+ ^ I [ (-2C^Cj^ + (-4 + -S{S + 1)) Cl^ Po 
- \c\c, + (-g + ^S{S + 1)) c^cl - 

X In^ iCpas) + Ci In (Cpo;,) + . . . | + . . . , 



(5) 



where Pi is the (i + l)-loop coefficient of the QCD /? function (/3o = 11Ca/12 — Tirn//3, . . .) 
and aj parameterizes the i-loop correction to the Coulomb potential (ai = 31Ca/9 — 
2QTpni/9, . . . [21]). For the processes under consideration, the total spin of the quark- 
antiquark pair is = 1. Nevertheless, we retain the full S dependence, so that our result 
is also applicable to processes with = 0, such as the decay rjb — > 77 or the production 
process 77 — at a future high-energy photon collider, which is dominated by the S 
wave [22]. The corrections through 0{a'^) have been derived in Ref. [9] for arbitrary n. 
The 0{a^\n^ as) correction has been obtained in Ref. [11]. In this Letter, we present the 
^(a^lnQ;^) correction by specifying the missing coefficient Ci in Eq. (5). 

The origin of the logarithmic corrections is the presence of several scales in the thresh- 
old problem. A logarithmic integral between different scales yields a term proportional 
to Inf. which becomes Ina^ for bound states that are approximately Coulombic, so that 
V oc ag. In effective-theory calculations, the scale defining the upper (lower) limit of a log- 
arithmic integral is set to infinity (zero), which induces a UV (IR) divergence. Thus, the 
logarithmic corrections can be identified with the effective-theory singularities, which dra- 
matically simplifies the calculation. Our analysis proceeds along the lines of Ref. [23] (see 
also Ref. [24]), where similar corrections have been considered for the QED bound-state 
of positronium. In the calculation, we employ the N^LO effective Hamiltonian derived 
in Ref. [7] and take into account the retardation effects due to the chromoelectric dipole 
interaction of the heavy-quark-antiquark pair with the dynamical ultrasoft gluons studied 
in Refs. [7,25]. Our result reads 



Ci 



-3 + 



CaCf + 



xS{S+l))cl 



Po + 




-CaCj.+ (^--+-S{S + 1) 



1 / 59 \ 
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xS{S+l)]C''}, + {-^+2\n2 + {l -In 2) S{S + 1) 



15 



X Cj,TF +^CACFTFni +(^- ^S{S + 1)) C|7>ni. 



(6) 



For the analysis of T sum rules and tt threshold production, one needs the extension of 
this result to arbitrary n, which, leaving aside the trivial n dependence of lip, 
given by 



c 



IS 



C„ = Ci + \-2CaCf + (^-4 + -S{S + 1) j C2 J /3o 
X ( Inn + ^i{n + 1) - 2n^2H ~^ + ^^ + y + ;^ 



+ 



]CICf + 



41 7 
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S{S + 1) CaC|, + 3C|, 



Inn 



(7) 



where \l'„(a;) = d"'lnr(a;)/(ia;", r(a;) is Euler's T function, and 7^ = 0.577216 ... is Euler's 
constant. 

The structure of the IR singularities in the Wilson coefficients can be read off from 
the UV-singular part of the effective-theory result for the wave function at the origin. In 
this way, we find the /j, dependence of Cy{iJ,) to be 



In 4 + ^*1 

mi TT 



(2) 



+ . . . 



(8) 



where the two- [20] and three-loop anomalous dimensions of the nonrelativistic vector 
current are 



7f = -IcaCf + (-1 + ls{S + 1)) Cl, 



(9) 



7f = + (-1 + 1)) c.cl - , 



IV 



^CaCf + ( 2 - ^S{S + 1)]C'f 
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+31n2 - ls{S+ 1)) Cl + {-l + ^ln2+ -^ln2) 

\ 4q /A 5 

X Sis + 1) jC7|r^ + -CaCfTfUi + (g - ^ 

x5(5+l))c|Tfn,, (10) 

with S — 1. We retain the full 5" dependence in Eqs. (9) and (10) because, for 5" = 
0, they give the anomalous dimension of the nonrelativistic pseudoscalar current ip^x^ 
which is relevant for two-photon processes. Note that the S* = result for Ci and 7^*^^^ 
corresponds to the d-dimensional spinor algebra in a regularization scheme adopted in [26] 
for the calculation of the two-loop hard corrections to the two-photon heavy quarkonium 
production and annihilation. 

Let us now explore the numerical significance of our results. Putting everything to- 
gether and substituting constants by their approximate numerical values, Eq. (2) becomes 

Ti fti (1 - 1.70as{mb) - 7 .98 a'^,{mb) + . . .) 
X {1 -0.30«, - 5.19a^ln«3 + 17.16 

- 14.38 In^ a, - 0.165 In a, + . . .) . (11) 

Note that the NNLO contribution due to the second term in the square brackets of Eq. (2) 
is included in the third factor on the right-hand side of Eq. (11), as the corresponding 
appearance of is of nonrelativistic origin and normalized at the soft scale. Evaluating 
this using as{Mz) — 0.1185 and rub = 5.3 GeV [8], we obtain 

Ti ^ rl'0(i -0.449nlo +1.771nnlo -0.766n3lo' +■ ■ ■)■ (12) 

where only the logarithmic N^LO terms are retained, which is indicated by the prime on 
the subscript. A similar analysis of Eq. (3) with rrit = 174.3 GeV yields 

Ri « R\^{1 -0.244nlo +0.438nnlo -0.196n3lo' + - ■ ■)■ (13) 

Without the 0(af Ina^) term, the N^LO contributions in Eqs. (12) and (13) read -0.560 
and —0.148, respectively. We learn the following: (i) while the coefficients in Eq. (11) 
and the analogous series for Ri sharply increase in magnitude as we pass from NLO to 
NNLO, this disquieting trend discontinues as we move on to N^LO; (ii) the coefficients 
of the known NNLO and N^LO terms are typically of order 10. These observations 
suggest that the magnitude of the coefficient of the missing non-logarithmic 0{al) term 
is unlikely to exceed this characteristic benchmark by far. This term would then be 
expected to yield corrections of order 25% and 3% to Fi and Ri, respectively. This 
provides us with an estimate of the residual uncertainty of our approximation as far as 
the perturbative QCD corrections are concerned. Moreover, the absence of a rapid growth 
of the coefficients along with the alternating-sign character of the series suggest that the 
higher-order corrections are hkely to be below these estimates. These observations can be 
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substantiated by investigating the scale dependence of Fi and Ri. In fact, the shifts in 
these quantities due to a variation of Hs by a factor of 2 are reduced from 50% and 13% to 
19% and 9%, respectively, as we pass from NNLO to N^LO. The latter values are in the 
same ballpark as the theoretical uncertainties estimated above. This renders a reliable 
perturbative description of heavy- quarkonium production and annihilation feasible. 

To conclude, we computed the logarithmically enhanced N^LO corrections to the 
heavy-quarkonium production and annihilation rates. Our results provide a useful hint 
on the general structure of the high-order corrections and open a new perspective for 
the theoretical analysis. Together with the 0{almq) result for the heavy-quarkonium 
spectrum [7,8], they constitute central ingredients for the high-precision analysis of T 
sum rules and tt threshold production in e'^e~ and 77 scattering. Calculation of the 
remaining non- logarithmic N^LO term appears to be mandatory for reducing the theoret- 
ical uncertainty further. Another challenging problem is to complete the resummation of 
the next-to-next-to-leading logarithms [13]. Our derivation of the three-loop anomalous 
dimension of the nonrelativistic vector current marks a major step in this direction. 

A.A.P. acknowledges discussions with A. Hoang. This work was supported in part 
by DFG Grant No. KN 365/1-1 and BMBF Grant No. 05 HTlGUA/4. The work of 
V.A.S. was supported in part by RFBR Project No. 01-02-16171, Volkswagen Foundation 
Contract No. 1/77788, and INTAS Grant No. 00-00313. 

Note added: 

In a recent paper [27] the infrared structure of the three-loop corrections to the heavy 
quark pair production current was analyzed in the effective theory framework. On the 
basis of this investigation the next-to-next-to leading logarithmic corrections to the heavy 
quarkonium threshold production were partially resummed. By expanding this result in 
as the 0{a^\nas) correction to the heavy-quarkonium production and annihilation rates 
can be obtained. The result announced in [27] for the S — 1 spin-one process agrees with 
our result for Ci, Eq. (6). 

References 

[1] T. Appelquist and H.D. Politzer, Phys. Rev. Lett. 34, 43 (1975). 

[2] V.A. Novikov, L.B. Okun, M.A. Shifman, A.l. Vainshtein, M.B. Voloshin, and V.l. 
Zakharov, Phys. Rev. Lett. 38, 626 (1977); 38, 791(E) (1977); Phys. Rep. C 41, 1 
(1978). 

[3] V.S. Fadin and V.A. Khoze, Pis'ma Zh. Eksp. Teor. Fiz. 46, 417 (1987) [JETP Lett. 
46, 525 (1987)]. 

[4] M.B. Voloshin, Nucl. Phys. B154, 365 (1979); Yad. Fiz. 36, 247 (1982) [Sov. J. Nucl. 
Phys. 36, 143 (1982)]; H. Leutwyler, Phys. Lett. 98B, 447 (1981). 

[5] I.Y. Bigi, M.A. Shifman, and N. Uraltsev. Annu. Rev. Nucl. Part. Sci. 47, 591 (1997). 



8 



[6] M.J. Strassler and M.E. Peskin, Phys. Rev. D 43, 1500 (1991). 

[7] B.A. Kniehl, A. A. Penin, V.A. Smirnov, and M. Steinhauser, Phys. Rev. D 65, 
091503(R) (2002); Nucl. Phys. B635, 357 (2002). 

[8] A. A. Penin and M. Steinhauser, Phys. Lett. B 538, 335 (2002). 

[9] J.H. Kiihn, A. A. Penin, and A.A. Pivovarov, Nucl. Phys. B534, 356 (1998); A. A. 
Penin and A.A. Pivovarov, Phys. Lett. B 435, 413 (1998); Nucl. Phys. B549, 217 
(1999); B550, 375 (1999); K. Melnikov and A. Yelkhovsky Phys. Rev. D 59, 114009 
(1999). 

[10] A.H. Hoang, M. Beneke, K. Melnikov, T. Nagano, A. Ota, A.A. Penin, A.A. Pivo- 
varov, A. Signer, V.A. Smirnov, Y. Sumino, T. Teubner, O. Yakovlev, and A. 
Yelkhovsky, Eur. Phys. J.direct C 3, 1 (2000). 

[11] B.A. Kniehl and A.A. Penin, Nucl. Phys. B577, 197 (2000); A.V. Manohar and I.W. 
Stewart, Phys. Rev. D 63, 054004 (2001). 

[12] M.E. Luke, A.V. Manohar, and LZ. Rothstein, Phys. Rev. D 61, 074025 (2000). 

[13] A.H. Hoang, A.V. Manohar, LW. Stewart, and T. Teubner, Phys. Rev. Lett. 86, 
1951 (2001); Phys. Rev. D 65, 014014 (2002). 

[14] A. Pineda, Phys. Rev. D 65, 074007 (2002); Phys. Rev. D 66, 054022 (2002). 

[15] A. Pineda and J. Soto, Phys. Lett. B 420, 391 (1998); Phys. Rev. D 59, 016005 
(1999); A. Czarnecki, K. Melnikov, and A. Yelkhovsky Phys. Rev. A 59, 4316 (1999); 
M. Beneke, A. Signer, and V.A. Smirnov, Phys. Lett. B 454 (1999) 137. 

[16] W.E. Caswell and G.P. Lepage, Phys. Lett. B 167, 437 (1986); G.T. Bodwin, E. 
Braaten, and G.P. Lepage, Phys. Rev. D 51, 1125 (1995); 55, 5853(E) (1997). 

[17] A. Pineda and J. Soto, Nucl. Phys. B (Proc. Suppl.) 64, 428 (1998); N. Brambilla, 
A. Pineda, J. Soto, and A. Vairo, Nucl. Phys. B566, 275 (2000). 

[18] M. Beneke and V.A. Smirnov, Nucl. Phys. B522, 321 (1998); V.A. Smirnov, Ap- 
plied Asymptotic Expansions in Momenta and Masses, ( Springer- Verlag, Heidelberg, 
2001). 

[19] S. Titard and F.J. Yndurain, Phys. Rev. D 51, 6348 (1995); A. Pineda, Nucl. Phys. 
B494, 213 (1997). 

[20] A. Czarnecki and K. Melnikov, Phys. Rev. Lett. 80, 2531 (1998); M. Beneke, A. 
Signer, and V.A. Smirnov, ibid. 80, 2535 (1998). 

[21] Y. Schroder, Phys. Lett. B 447, 321 (1999) and references cited therein. 



9 



[22] A.A. Penin and A.A. Pivovarov, Nucl. Phys. B550, 375 (1999); Yad. Fiz. 64, 323 
(2001) [Phys. Atom. Nucl. 64, 275 (2001)]; A. Czarnecki and K. Melnikov, Phys. 
Rev. D65, 051501 (2002). 

[23] B.A. Kniehl and A.A. Penin, Phys. Rev. Lett. 85, 1210 (2000); 85, 3065(E) (2000); 
85, 5094 (2000). 

[24] R.J. Hill and G.P. Lepage, Phys. Rev. D 62, 111301 (2000); K. Melnikov and A. 
Yelkhovsky, ibid. 62, 116003 (2000); Phys. Rev. Lett. 86, 1498 (2001); A.A. Penin, 
Nucl. Phys. B (Proc. Suppl.) 96, 418 (2001); R.J. Hill, Phys. Rev. Lett. 86, 3280 
(2001). 

[25] B.A. Kniehl and A.A. Penin, Nucl. Phys. B563, 200 (1999). 
[26] A. Czarnecki and K. Melnikov, Phys. Lett. B 519, 212 (2001). 
[27] A. Hoang, Report No. MPP-2003-38 and hep-ph/0307376. 



10 



